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Abstract 

Using a sample of 3.3 x 10^ T{4S) BB events collected with the CLEO II 
detector at the Cornell Electron Storage Ring (CESR), we measure B{B 
piv), \Vub\, and the partial rate (AF) in three bins of = {pB—Ppf'- We find 
B{B^ p-i+v) = (2.69 ± 0.4lt° j5 ± 0.50) x lO'^, |Kfe| = (3.23 ± 0.24l°:23 ± 
0.58) X 10-3, AF(0 < g2 < 7 GeVVc^) = (7.6 ± S.Ot?;^ ± 3.0) x lO'^ ns-\ 
AF(7 < g2 < 14 GeVVc^) = (4.8 ± 2.9t[}:^ ± 0.7) x lO^^ ns-\ and AF(14 < 
< 21 GeVVc^) = (7.1 ± 2.ltil ± 0.6) x 10"^ ^^-i ^ ^ jjere, i = e ox p, but 
not both, and the quoted errors are statistical, systematic, and theoretical. 
The method is sensitive primarily to i? — > p^v decays with leptons in the 
energy range above 2.3 GeV. Averaging with the previously published CLEO 
results for B piv, we obtain B{B^ p-i+v) = (2.57 ± 0.29t|];^^ ± 0.41) x 
10-4 and \Vub\ = (3.25 ± O.Utg:^! ± 0.55) x lO^^. 
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I. INTRODUCTION 



Exclusive semileptonic b — > uiu decays are an active area of experimental and theoretical 
study p|-p3|. These rare processes can be used to extract the magnitude of Vub, one of the small- 
est and least well known elements of the Cabibbo-Kobayashi-Maskawa (CKM) quark-mixing ma- 
trix [24|. Because ~ 0.08, the branching fractions for exclusive b — > uii^ processes are 



small, of order 10~^, and they have only recently become experimentally accessible. The present 
analysis confirms the initial CLEG observation [|^ of — > pii/, improves the precision on the 
branching fraction and \Vub\, and provides the first information on the dependence for this 
decay. 

Extracting \ Vub\ from a measured decay rate requires significant theoretical input because the 
matrix elements for such processes involve complex strong-interaction dynamics. Although the 
underlying b — > uii' decay is a relatively simple weak process, it is difficult to calculate the strong- 
interaction effects involved in the transition from the heavy B meson to the light daughter meson. 
Because of these theoretical uncertainties, even a perfectly measured B — > pii' branching fraction 
would not at present lead to a precise value of \ Vub\- 

The dynamics in B ^ piu decay are in contrast with b — > civ decays, such as B ^ D*iu, 
where a heavy quark is present both in initial and final states. In this case, techniques based on 
Heavy Quark Effective Theory (HQET) can be used to calculate the decay amplitude with good 
precision, particularly for the kinematic configuration in which the charm hadron has zero recoil 
velocity. The zero-recoil point in i? — > plv cannot be treated with similar techniques, however, 
because the daughter u quark is not heavy compared to the scale of hadronic energy transfers. 
Nevertheless, substantial progress has been made using a variety of theoretical methods, including 
quark models |6-12], lattice QCD [^^-|15|, QCD sum rules |p!^18|, and models relating form factors 



measured in D ^ K*£u decay [19| to those in B ^ p£u decay. 

Experimentally, the main difficulty in observing signals from b — > uii' processes is the very 
large background due to 6 ^ cii'. Because a significant fraction of i? ^ piv events have lepton 
energy beyond the endpoint for b civ decay, lepton-energy requirements provide a powerful tool 
for background suppression. However, extrapolation of the decay rate measured in this portion of 
phase space to the full rate again requires the use of theoretical models, and it introduces model 
dependence beyond that associated with simply extracting the value of \Vub\ from the branching 
fraction. In this study, we begin to explore the decay dynamics oi B ^ plv by measuring the 
distribution of q^, the square of the mass of the virtual W . The distribution of q^ is reflected in 
the p momentum spectrum. Eventually, studies of the cf' distribution, as well as of the angular 
distributions of the decay products, should reduce the model dependence on \Vub\ by constraining 
theoretical models for the decay form factors. 

CLEG has previously observed the decays 5° p~i^i^, B+ p^i+v, B° tt'I+u, 
and B^ — > Tr^i'^v by measuring both the missing energy and momentum in an event to infer the 
momentum of the neutrino. Using 2.84 x 10^ T(45') — > BB events, this study obtained 

B{B^ p-e+u) = (2.5 ± 0.4l°:^ ± 0.5) x 10-^ 
^ ir'i+i^) = (1.8 ± 0.4 ± 0.3 ± 0.2) x 10"^ 

\Vub\ = (3.3 ± 0.2^0:1 ± 0.7) X 10-3, (1) 

where the errors are statistical, systematic, and theoretical. A large contribution to the B piv 
systematic error is associated with a possible nonresonant B vrvr^z/ rate. In the analysis described 
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below this uncertainty is reduced. 

We report on a measurement of B{B ^ piv) and \Vub\ as well as the first measurement of the 
partial rate for B — > piu in bins of q'^. We study five signal modes: p~i~^iy', B~^ p^i~^u, 

B~^ — s- ojI'^u^ B^ — > 7r~£"'"z^, and B^ — > tt^I'^u. Our method is sensitive primarily to leptons with 
energies above the h — > civ lepton-energy endpoint {E^^ > 2.3 GeV). The resulting event sample is 
essentially statistically independent of that from the previous analysis since it contains much larger 
signal and background yields, even though the data samples are similar. 

In this paper, we use the notation b u{c)iv to denote inclusive B X^(^^-^lv decay, where 
Xu is a hadronic system with no charm quarks, and Xc is a hadronic system with a charm quark. 
When discussing backgrounds, we will often refer to 6 ^ ulv inclusive decays excluding one or more 
signal modes. In this case, we will use the notation h — > ulv and explicitly list the modes not to be 
included. When no charge state is given, B — > piv refers generically to 5 — > p^iv, and B p^iu 
decays. When specifying a particular decay mode, we implicitly include its charge-conjugate decay. 

We begin in Sec. ^ by introducing the phenomenology of -B ^ piu decay. The full differential 
decay rate is expressed in terms of helicity amplitudes, which depend in turn on form factors, 
Lorentz- invariant functions of q'^ that parametrize the hadronic current. We discuss five form- 
factor models that are used to obtain ll^^l and to extrapolate our yield measured at high lepton 
energy to the full phase space available in i? — > piv decay. 

Section [D describes the data sample and the requirements used to distinguish signal events 
from backgrounds, which are due primarily to continuum events (e"'"e~ — > qq, qq = uu, dd, ss, and 
cc), b — > cii', and b uiv (other than B plv or i? — > u^lv) events. 

Section IV describes the binned maximum-likelihood fit used to extract the B — > p(.v signal. 
Although we must rely on Monte Carlo calculations to model the shapes of distributions for the 
b — c{u)iv backgrounds, the normalization and lepton-energy spectrum of each background com- 
ponent are determined by the data. In addition, the assumed background shapes are extensively 
tested using sideband regions where the contribution from signal events is small. The continuum 
background contribution is measured directly using data. Section ^ presents the yields extracted 
from the fit and the model-dependent detection efficiencies used to obtain values for B{B plu) 
and iKbl- 

Section VI describes the contributions to our systematic error. Apart from model dependence, 
these are due primarily to uncertainties in the b uiu backgrounds and to uncertainties in the 
detector simulation. In the high lepton-energy region, where the sensitivity to the signal is greatest, 
the background from other b — uiiy processes is comparable to that from b — > ciu decays, but its 
kinematic properties are less well understood. 

Section VII presents our measured values of B{B ^ plu) and \ Vub\i as well as the average with 
the previous CLEO result. This average takes into account the correlation in the systematic errors 
as well as the variation in results for different form-factor models. In Sec. VIII we present our 
conclusions and an outlook for future measurements. 



II. SEMILEPTONIC DECAY KINEMATICS 

The matrix element for the semileptonic decay P{Qq) — > V{q'q)i~i' of a pseudoscalar meson, 
P(Qq), to a vector meson, V{q'q), can be written 



M{P ^ Viu) = -i^Vg^QL^J^, (2) 
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where Vq/Q is the CKM matrix element for the Q ^ q' transition L'^ is the leptonic current 



(3) 



and the hadronic current 



= {V{pv,e)\{V - A)^\P{pp)) = - e* (mp + mv)Ai{q^) 

+ {pv +pp)^{e* ■ q) 



[mp + my) 



+ q^.{e*-q)'^(.A,{q^)-Ao{q^)) 



+ "^i^fiupa^ PyPp 



mp + mv ' 



(4) 



is written in terms of four form factors, Ai{q^), A2{q^), V{q^), and AQ{q^), where ^o(O) = ^3(0) 
and 



A,{q^ 



mp + mv ^ , 2n mp - my 
2my ^ 2my 



A2{q^). 



(5) 



Here mp, my, pp, and py are the pseudoscalar and vector meson masses and four-momenta, 
q = Pp — Pv, and e is the vector- meson polarization four- vector. Terms in proportional to q^ 
vanish in the limit of massless leptons, so that the decay B piu depends effectively on only three 
form factors (yli(g^), A2{q'^), and V{q'^)) for electrons or muons. 

The amplitudes for the vector meson to have helicity +1, —1, or 0, denoted by H+{q'^), H^{q'^), 
and Ho{q'^), can be expressed in terms of the form factors: 

H^iq') = {mp + my)A,iq') T ^^^^^V{q'), 



1 



2myyjq- 



mp + my 
:[(mp - my - q^){mp + my)Ai{q^) 



^ mp\PV? 
mp + my 



A2{q^)], 



(6) 



where py is the vector meson three-momentum in the 5-meson rest frame, which is a function of 
^2 



iPvl 



\ 



mp + my 
2m p 



mr 



As q'^ decreases from q'^ 



{mp — myY to 



0, |py| increases from |py| 



(7) 



to 



(mp — my) /2m^p. At large (f' , the A\{(f') form factor therefore dominates all three of 



iPvlma 

the helicity amplitudes. The full differential decay rate for P Viu {V — > P1P2) is given by |25 

mi, 



dq'^ d cos 9 id cos By dx 8{4tt)'^ 

X {l-7]cos9efsm'^9y\H+{q'^)f 
+ (1 + r/ cos 0^)2 sin^ 9y\H_{q'^)\'^ 
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FIG. 1. Model predictions for the (a) ^i(g^), (b) ^2(9^), and (c) V{q'^) form factors. The 
Ai{q'^) form factor determines the rate at high q^. The V{q'^) form factor has a large effect on 
the lepton-energy spectrum. Maximum daughter-hadron recoil occurs at the minimum value of q"^, 
which is close to zero. 



+ 4sin2 0^cos2 ev\Ho{q'^)f 

— 4r] sin 0^(1 — i] cos Of) sin Oy cos Oy cos xH+{q'^)Ho{q'^) 
+ At] sin 0^(1 + 7] cos 0£) sin Oy cos 6y cos xH-{q'^)HQ{q'^) 

- 2sin2 6'^sin2 Oy cos2xH+{q^)H^{q^) 



(8) 



where q'^ = {pp — py)"^, Oi is the polar angle of the lepton in the W rest frame with respect to 
the W flight direction in the B rest frame, and 9y is the polar angle of one of the pseudoscalar 
daughters in the rest frame of the vector meson with respect to the vector-meson flight direction 
in the B rest frame, x is the angle between the decay planes of the W and the vector meson. The 
factor r] is equal to +1 (—1) when the quark Q has charge —1/3 (+2/3). 

For the extraction oiB{B — > piv) and AF in bins of q^ from the data, we use Eq. ^ to extrapolate 
our measurements based primarily on the high lepton energy part of the decay phase space to 
the entire phase space. The form factors ^i(g^), A2(g^), and T^(g^) are evaluated using various 
models, whose predictions are shown in Fig. ||. Models additionally provide the normalization 
Fthy = F(i? ^ /j£z^)/|T4bp, which allows us to extract \Vub\ from the measured decay rate. We use 
a set of form-factor models representative of current theoretical work, including results from two 
quark-model approaches (ISGW2 |] and Melikhov/Beyer 0), lattice QCD (UKQCD light- 
cone sum rules (LCSR [^), and a method incorporating form-factor measurements in D K*iv 
decay (Wise/Ligeti+E791 In all cases we have used form-factor parametrizations as described 

in the above references. A theoretical error is added to the branching fraction and AF measurements 
based on one-half of the full spread in results among these models. The \Vub\ measurement has an 
additional theoretical uncertainty due to the determination of Fthy 

The Dalitz plot for B — > plv decay predicted by the ISGW2 model is shown in Fig. ^. For 
a given value of g^, cosO^ varies from —1 to +1 across the allowed range of lepton energies. The 
Dalitz plot shows that positive values of cos^£ are favored over negative values; this effect, which 
produces a hard lepton-energy spectrum, is a consequence of the V — A coupling, and is discussed 
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FIG. 2. A simulated Dalitz plot, vs. E'^, for B — > piv decay. The lepton energy, E^, is 
computed in the T(4S') — > BB rest frame (lab frame), which differs slightly from the B rest frame. 
This analysis is most sensitive to i? — > p^v events in the lepton-energy region to the right of 
the vertical line, > 2.3 GeV. The horizontal lines define the bins in the AT measurement: 
< g2 < 7 GeVVc^, 7 < g2 < 14 GeVVc^ and U < q^ < 21 GeVVc^- At high our 
lepton-energy requirement retains a much larger fraction of the phase space than it does at low q^ . 



further below. We are most sensitive to S — > plv decay in the endpoint region to the right of the 
vertical line, E(_ > 2.3 GeV. Although a significant fraction of the total predicted rate for B — > plu 
lies at high lepton energy. Table || shows that theoretical approaches differ, giving results between 
24% and 35% for this fraction. 

Table | also shows the quantity Fthy predicted by each form- factor model. Even with a perfectly 
measured B piu branching fraction, \Vub\ would be uncertain due to the spread in T^hy values 
among form-factor models. We can summarize the model predictions as Tthy = (16.8 zt 2.6) ps^^, 
a 15% spread. As explained in Section VII, we assign a 30% error on Tthy 

Figure |^ shows the distribution predicted by our implementation of these models. Integrated 
over all lepton energies, there is a substantial difference among the models at low (large daughter- 
hadron recoil). At high lepton energy, however, the shapes of the q"^ distributions predicted by each 
of the form-factor models are very similar. To distinguish among form-factor models on the basis 
of the distribution, the lepton-energy region below 2.0 GeV must be probed. 

Table || shows how the efficiency of the lepton-energy requirement varies with for each form- 
factor model. Our lepton-energy requirement retains a much larger fraction of the phase space at 
high q'^ than it does at low g^. Thus, at high g^, the efficiency is quite high and a small variation 
is seen among models. At low g^, on the other hand, the efficiency is relatively low with a large 
variation among the form-factor models. 

Integrating Eq. S over the angular variables, we obtain 



\H^{q')\' + \H^{q')U\Hoiq')\' 



(9) 
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FIG. 3. Theoretical distributions for the B piu form-factor models that we have consid- 
ered: (a) over the full lepton-energy range and (b) for En > 2.3 GeV. The lepton energy, Ei, is 
computed in the T(45) — > BB rest frame (lab frame), which differs slightly from the B rest frame. 
In each plot, models have been normalized to equal areas. For high lepton energies, Ei > 2.3 GeV, 
all of the models predict very similar distributions. 



FF model Fthy (ps~^) TjEj > 2.3 GeV)/r (%) r(2.0 < E^ < 2.3 GeV)/r (%) 



ISGW2 14.2 35 33 

LCSR 16.9 24 28 

UKQCD 16.5 27 30 

Wise/Ligeti-hE791 19.4 31 34 

Beyer/Melikhov 16.0 27 30 



TABLE I. Fthy and lepton-energy distribution predictions from each form-factor (FF) model. 
The lepton energy, Ei, is computed in the T(45') BB rest frame (lab frame), which differs 
slightly from the B rest frame. For the branching fraction measurement, the primary source 
of model dependence is in the ratio T(E£ > 2.3 GeV)/r. For the \Vub\ measurement, there is 
significant additional model dependence arising from the quantity Fthy 
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FIG. 4. dV /dq^ distributions for each of the three terms in Eq. |9|: (a) the terms proportional 
to l-ff-P and l-f^+P and (b) the |i^oP term. 



Figure |^ shows the dT/dq"^ distribution for each term in Eq. |^. In each form- factor model, the 
term proportional to contributes the largest fraction of the total rate. As a consequence of 

the V — A couplings of the W boson the vector meson is more likely to have helicity —1 than +1. 
For B — > piu decay, this asymmetry is quite large, except at very high q^ (low daughter-hadron 
recoil momentum) where the Ai form factor dominates each of the helicity amplitudes and the p is 
nearly unpolarized, and at small values of q^, where only the zero- helicity component contributes 
significantly. Because the contribution is quite small, the theoretical error in the B{B piv) 

and AF measurements is due primarily to the uncertainty in the relative size of H-{(f') and HQ{q^). 

Table ^ shows the predictions for F/|14b|^ in bins of q^ (AFthy) for the form-factor models. 
These predictions can be used to extract from a measurement of the q^ distribution. 



III. DATA SAMPLE AND EVENT SELECTION 

The data used in this analysis were collected using the CLEO II detector [27| located at the 



Cornell Electron Storage Ring (CESR), operating near the T(4S) resonance. We have analysed a 
3.1 fb~^ sample taken on the resonance, corresponding to approximately 3.3 x 10® BB pairs. Addi- 
tionally, we examine a 1.7 fb~^ sample taken at a center-of-mass energy 60 MeV below the T(4S') 
mass. These off-resonance data, which have an energy below the threshold for BB production, are 
used to measure the continuum background. 

The CLEO II detector is designed to provide excellent charged- and neutral-particle recon- 
struction efficiency and resolution. Three cylindrical tracking chambers are surrounded by a time- 
of-flight (TOF) system, a Csl calorimeter, and a muon tracking system. The nearly complete 
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FF model 


Ar{Ei > 2.3 GeV)/Ar 

(%) 

OKq"^ <7 GeVVc^ 


Ar(£;^ > 2.3 GeV)/Ar 

(%) 

7 < g2 <- 14 GeVVc^ 


Ar{Ei > 2.3 GeV)/Ar 

(%) 

14 < g2 < 21 GeVVc^ 


ISGW2 


24 


35 


43 


LCSR 


12 


25 


38 


UKQCD 


15 


28 


39 


Wise/Ligeti+E791 


23 


29 


37 


Beyer/Melikhov 


14 


27 


38 



TABLE II. The fraction of events predicted to have > 2.3 GeV in bins of for each 
form-factor (FF) model. The spread among these results is the largest component of the theoretical 
error in the AF measurement. 



FF model 


Afthy (pS-1) 

< g2 < 7 GeV7c^ 


Afthy (pS-1) 

7 < g2 < 14 GeVVc^ 


Afthy (PS-^) 

14 < g2 < 21 GeVVc^ 


ISGW2 


2.8 


6.1 


5.3 


LCSR 


5.4 


6.6 


4.9 


UKQCD 


4.5 


6.7 


5.3 


Wise/Ligeti+E791 


3.5 


8.4 


7.5 


Beyer /Melikhov 


4.3 


6.4 


5.3 



TABLE III. Partial rate (Afthy) predictions in bins from each form-factor (FF) model. 
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solid-angle coverage of both the tracking system and calorimeter is exploited in this analysis to 
obtain information on the momentum of the neutrino. 

To be considered as a hadron or lepton candidate, charged tracks must satisfy several require- 
ments. The impact parameter of the track along and transverse to the beam direction must be less 
than 5 cm and 2 mm, respectively. The spread in the impact parameter along the beam direction 
is dominated by the beam width in that direction and is approximately 1.5 cm. Transverse to 
the beam direction, the resolution on the impact parameter is approximately 0.3 mm. The r.m.s. 
residual for the hits on the track must be less than 1 mm. Finally, we require that at least 15 out 
of the 67 tracking layers be used in the track fit. 

Electron and muon identification requirements are chosen to achieve high efficiency and a low 
hadronic fake rate. We require | cos 9\ < 0.85, where 6 is the polar angle of the track's momentum 
vector with respect to the beam axis. Electrons are identified primarily using the ratio of the 
calorimeter energy to track momentum (E/p) and specific ionization (dE/dx) information. Muon 
candidates are found by extrapolating tracks from the central detector into the muon counters; such 
candidates are required to penetrate at least seven interaction lengths of iron. For electrons (muons) 
that satisfy the tracking requirements described above as well as the lepton-energy requirements of 
our analysis, the efficiency of these lepton-identification requirements is approximately 94% (82%). 
The fraction of hadrons passing electron (muon) identification criteria is determined from data and 
is approximately 0.1% (0.6%). 

Photon candidates are associated with Csl calorimeter clusters that are not matched to any 
charged track. Individual photons are required to have energy, Ej, greater than 30 MeV. To be 
considered in a vr^ candidate, two photons must have a combined invariant mass, M^^, within 2a 
of the nominal vr'' mass. Their combined energy, E^^, is required to be greater than 325 MeV. In 
addition, we require that at least one photon satisfy |cos6'-y| < 0.71, where 9^ is the polar angle 
with respect to the beam axis. Each photon candidate is included in at most one 7r° candidate. 

To search for B piu, B coii', and B — > Triu events, electron or muon candidates are com- 
bined with TT^TT^, TT+TT", 7r"'"7r~7r'^, TT^, Or TT^ Candidates. Combinatoric backgrounds are suppressed 
by requiring candidates for these hadronic systems to have a momentum greater than 300 MeV/c; 
these backgrounds are further suppressed in the tt'^tt^tt^ mode by making a cut on the io 7r"'"7r~7r'^ 
Dalitz amplitude. We require that the it'^,it~, and 7r° be configured such that the decay probability 
density {dV oc |p^+ x p^^- P where Pt^+ and p^^- are evaluated in the co rest frame) is at least 12% 
of its maximum value. 

We define the missing momentum in the event to be 

Pmiss = -Spi, (10) 

where the sum is over reconstructed charged tracks and photon candidates in the event. If the 
event contains no undetected particles other than the u from the B p£u decay, Pmiss ~ Pu- The 

resolution on Pmiss is determined by undetected particles, such as mesons, and our experimental 
rcsohition. As discussed in the following section, Pmiss is used to calculate one of our three fit 
variables. 

Background events come from continuum, b — > ctv, b uiv (other than B — > piu or B ^ Loiv), 
and events with hadrons misidentified as leptons. The b civ background is dominant at lower 
lepton energy. To distinguish signal events from background most effectively, we consider three 
lepton-energy ranges: 

• Ei > 2.3 GeV (denoted as HILEP in the remainder of this paper). 
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• 2.0 <Eii< 2.3 GeV (LOLEP), 

• 1.7 <Ei< 2.0 GeV, 

where the lepton energy, E^, is computed in the T(4S') — > BB rest frame (lab frame). The S-meson 
rest frame is moving with a small velocity, (3 ~ 0.065, relative to the T(45') frame. Sensitivity to 
B — > plv decay is best in the HILEP lepton-energy bin where the h — > civ background is quite 
small. Here, the backgrounds are primarily continuum and h — > ulv events. In LOLEP, the b — > civ 
background dominates, but the h — > ulv contributions in this region are not negligible. Events with 
1.7 < El < 2.0 GeV are completely dominated by 6 — > civ decays and are included primarily to 
verify our understanding of this source of background. 

Because our best sensitivity to B ^ piv is in the HILEP region, suppression of continuum 
background is a key issue in the analysis. At the T(45) center-of-mass energy, T(4S') — > BB 
decays can be distinguished from continuum events, which have a three times larger cross section, 
on the basis of various quantities that describe the overall distribution of tracks and photons in 
the event. Such event shape, or event topology, variables exploit the fact that the B mesons are 
produced nearly at rest so their decay products are distributed roughly uniformly in solid angle. In 
contrast, continuum events yield a much more collimated or jet-like event topology. We use several 
event-shape requirements in selecting our sample: 

• The ratio of the second to zeroth Fox- Wolfram moments is required to be less than 0.3. 
This ratio of moments tends to be close to zero for spherical events and up to unity for jet-like 
events. 

• pi, uji, and 7ri candidates are required to have | cos^thrustl < 0.8, where ^thrust is defined to 
be the angle between the thrust axis of the candidate system and the thrust axis of the rest 
of the event. 

• We select events with | cos ^miss| < 0.9, where ^miss is the polar angle of the missing momentum 
(Pmiss) with respect to the beam axis (see Eq. |^). 

• We select events whose energy is evenly distributed around the momentum axis of the candi- 
date lepton using a Fisher discriminant ||2^. The input variables to the Fisher discriminant 
are the track and Csl cluster energies in nine cones of equal solid angle around the lepton- 
momentum axis. Energy in BB events tends to be more evenly distributed in these cones 
than in jet-like continuum events. 

To suppress combinatoric backgrounds from b civ and from b — > uiv sources other than 
B piv and B — > ujiv, we require that the event kinematics be consistent with these signal 
modes. Using the constraints Eb = EBea.m and = = {pb — Pp — Pif'^ we compute the angle 
between the B momentum direction and that of the reconstructed Y = (p, uj, it) + i system: 

_ 2EbEy -ml -ml 
cos 9 BY = ^1 — rr^ ■ (11) 

The distribution of \pb\ ■ cos 9 by is shown in Fig. |5|a for signal, continuum, and b — > civ events 
that satisfy all other analysis requirements. We require that \pb\ ■ \ cos9by\ ^ 385 MeV/c. For 
well-reconstructed signal events with a B momentum of less than 385 MeV/c, this cut is 100% 
efficient, but b uiv (non-signal modes), b civ, and continuum events often produce unphysical 
values of cos 9 by and can be rejected. At CESR, the average B meson momentum is approximately 
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FIG. 5. Signal and background distributions for the HILEP vrvr modes: (a 
for events that pass all other requirements including |M(7r7r) — M(p)\ < 0.15 GeV/c 
\AE\ < 0.5 GeV and (b) M(7r7r) for events that pass all requirements including lAi^j < 0.5 GeV. 
Shown are the Monte Carlo signal, continuum data, and Monte Carlo b — > cii' contributions. Each 



has been normalized using the results of the likelihood fit described in Sec. IV. The signal contri 



bution includes events from one signal mode that have been reconstructed in another signal mode. 
The \pb\ • I cos6by\ < 385 MeV/c requirement keeps almost all signal events, while rejecting much 
of the background. 
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320 MeV/c, with a standard deviation of 50 MeV/c due to the beam energy spread from the emission 
of synchrotron radiation, so this cut is nearly 100% efficient for well-reconstructed signal events. 

Finally, we compare the direction of the missing momentum (pmiss) with that of the neutrino 
momentum (pj^) inferred from the Pu = Pb — Py momentum. The latter is known to within an 
azimuthal ambiguity about the py direction because the magnitude, but not the direction, of the 
B meson momentum is known. We require that the minimum difference in the angle between these 
two directions be less than 0.6 radians. 

Because the p has a large width and we are reconstructing several different hadronic modes, 
events with multiple entries passing all of the above requirements are common. To avoid the 
statistical difficulties associated with this, we choose one combination per event after all other 
criteria are imposed, picking the combination with |py + PmissI closest to |p_b|. 

The large p width also leads to an important effect in signal events. The candidate vrvr system 
may not consist of the true daughter particles of the p. For example, it is possible for a S — > p^iu 
event to satisfy the analysis requirements for the B p'^iv channel. In addition, B ujii' events 
can feed into the B plv channels. We denote this contribution as crossfeed signal. Although the 
kinematic distributions of these events are somewhat different from correctly reconstructed signal 
events, the crossfeed contribution is produced by signal processes, and it is counted as such. 

Figure ^b shows the M(7r7r) distribution expected from signal, continuum, and b — > ciu events 
in the HILEP vrvr modes. Each contribution has been normalized using the results of the likelihood 
fit described in Sec. IV. 



IV. FIT METHOD 



We perform a binned maximum-likelihood fit using MINUIT |3C| in three lepton-energy ranges 
and five decay modes. In addition to lepton energy, we fit two kinematic variables, so we are in 
effect performing a three dimensional correlated fit. Our fit includes contributions from B — > piu, 
B Lviv, B Tilv, other h uiv modes, b — > cii', continuum events, as well as a contribution 
from hadrons faking leptons. In this section, we describe the fit variables and constraints, the 
modeling of signal and background components, and our method for extracting the rate as a 
function of q^. 

The three fit variables, Ei, M^ad) and AE, are constructed from the three decay products in 
semileptonic decay: the lepton, the final-state hadron, and the neutrino. The three lepton-energy 



bins were discussed in Sec. III. The invariant mass of the daughter hadron, Mhad, is reconstructed 
in the vr'^'vr" and vr^vr*^ channels for B plv and in the vr^vr^vr*^ channel for B ojlu. The bin 
size for the fit is 90 MeV/c^ for both M(7rvr) and M(7rvr7r). The momentum of the hadron is not 
used as a fit variable but is the basis for measuring the q'^ distribution after the fit is performed. 

The presence of a neutrino consistent with B — > piv decay is signaled by a peak in the distri- 
bution of 

AE = Ep + Ei + Ipmissi — -E'Seam (12) 

at AE = 0. Backgrounds from b civ and b uiv (except for B — > tt(.v) events tend to have 
AE < 0. On the other hand, B niv events have AE > when reconstructed in the vrvr or 
vr+vr^vr*^ modes, since extra particles beyond the actual B decay products are included in the decay 
hypothesis. The AE bin size for the fit is 200 MeV. 
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/S.E vs. M['KTi) distributions for signal and background contributions in the HILEP region are 
shown in Fig. ^. Shapes of distributions for signal, h — > civ, and h —>■ uiu are taken from Monte 
Carlo simulations using theoretical models of i?-meson decays and full detector simulations based 
on GEANT ||3l| , while the continuum background is measured with the off-resonance data. An ad- 
ditional background, not shown, arises from hadrons passing the lepton identification requirements 
in BB events. This contribution is determined using measured fake rates and data, as described 
below. The data or Monte Carlo samples that are used to estimate the non-continuum contribu- 
tions are sufficiently large that their statistical fluctuations are negligible compared with those of 
the data. We now discuss each of these fit components in more detail. 

For signal, the Monte Carlo sample provides the shapes of the three-dimensional distribution of 
El, Mhad) and AE. There is very little model dependence in the AE vs. Mhad distributions, but 
the distribution of signal events across the three lepton-energy bins varies significantly among form- 
factor models. The variation among the B — > plu form- factor models is a significant systematic 
error for our measurement of B{B — > piv) and \ Vub\- 

There are a large number of crossfeed signal events, as described in the previous section. The 
size of the crossfeed contribution relative to the direct signal is determined using Monte Carlo 
simulation. These events are misreconstructed signal decays and their kinematic distributions can 
be somewhat different from direct signal events. 

Isospin and quark-model relations are used to constrain the relative normalizations of B^ 
p^i^v, B^ p^t^v, and B^ — > ujt^v and, separately, those of B^ — > -K^i'^v and B^ — > ir^t^v. 
Isospin symmetry implies that 

r(5° ^ p-tu) = 2T{B+ p^i^iy), 

r(5° ^ 7:-£+iy) = 2r{B+ ^ 7r°^+i/), (13) 
while the u and p^ wave functions are expected to be very similar in the quark model, giving 

T{B+ p^l+u) = T{B+ Lol+v). (14) 



Possible isospin breaking effects are discussed in Refs. [32,3^. We assume that the T(4S') decays 
only to BB in equal proportions of B^B^ and B~^B~ mesons (/-| //oo = 1-0 =b 0.1). We also use 



the value tb+ /t^o = 1.04 it 0.04 [34| for the ratio of charged to neutral i?-meson lifetimes. 



The shape and normalization of the continuum background are measured using the off-resonance 
data sample. The statistical fluctuations in the off-resonance sample are accounted for by fitting the 
on- and off-resonance samples simultaneously, accounting for the difference in luminosity and cross 
section. For each bin in the fit, the mean number of continuum background events is determined 
so as to maximize the combined likelihood. 

The normalization of the h clu background is determined to a large extent from the data. 
This background arises primarily from two extensively studied decays, B Dlv and B D*iv, 
and the associated kinematic distributions are well known. As for the signal modes, the shapes of 
the AE vs. Mhad distributions due to 6 ^ ciu are obtained using the Monte Carlo calculations. 
With these constraints, there remain 15 6 ^ ciu normalization parameters, one for each of the five 
signal modes in each of three lepton-energy bins. 

In our fit, there are in fact only seven free parameters associated with the h — > cd-u background, 
defined in the following way. Five of the parameters, one for each signal mode, are scale factors 
that give the overall normalization of the h — > ciu background relative to that expected from the 
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FIG. 6. Signal and background AE vs. M(7r7r) distributions for the combined tt^vt" and 7r"'"7r~ 
modes in HILEP: (a) the direct signal contribution; background events from (b) the continuum; (c) 
b civ events; (d) h — > ulv events. The continuum contribution is modeled using the off-resonance 
data, while the other contributions are taken from Monte Carlo samples. Signal events are centered 
on /S.E = and are clustered around the p mass in M{ti'k). 
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Monte Carlo simulation. Two additional parameters describe the ratios of scale factors among the 
three lepton-energy bins. A common set of ratios is used for all five modes, so we have two rather 
than 2x5 parameters for the b ciu lepton energy spectrum. The reason for this choice is that 
the same set of background modes, dominated hy B ^ D*ii' and B Diu, contributes to each 
of the five signal modes. The relative contributions of the different b — > cii' modes are essentially 
independent of the b — > ulu mode we reconstruct. In Sec. 0, we show that our fit results for all 
seven of these b civ background parameters agree well with Monte Carlo predictions. 

The normalization for the contribution from B — > it(.v decay is also determined by the data. As 
for the B — > plv contribution, form-factor models are used to describe the kinematic distributions. 
A systematic error is assigned based on the variation observed in the B — > piu fit results for different 
B nil/ form-factor models. 

The ISGW2 Q predictions are used to model the distributions of 6 — > uiv sources other than 
B piu, B — > ujiv, and B — > ttIv for resonances up to the p(1450). The dominant resonances are 
r/, 6i, and oi. The uncertainty on the composition of this background is considered as a systematic 
error. To reduce this systematic uncertainty, we allow for a separate normalization for the b uiv 
component in each lepton-energy bin. 

The small contribution from fake leptons is found using the data. The shapes of the AE 
vs. M^tttt) distribution, along with that of E£, are found by combining p candidates with charged 
tracks not identified as leptons. These combinations are required to pass all analysis criteria 
except for the lepton-identification requirements. The normalization of the fake-lepton component 
is determined by measurements of the probability for hadrons to satisfy the lepton identification 
requirements using vr, K, and p samples from the data (see Sec. 

There are thus twelve free parameters in the fit: 

• B{B piu) (1 parameter). 

• B{B Triu) (1 parameter). 

• The yield of 6 — > uiu events in each lepton-energy bin (3 parameters). 

• The seven parameters that describe the normalization of the b civ contribution. 

We do not include a contribution from B Tnriv nonresonant events in our fit. Instead, we 



use the data to constrain the size of this contribution and assign a systematic error (see Sec. VI) 
The B{B — > p£v) result from the likelihood fit is used to compute 



V J-thyT'BO 

B is the result for the B^ p^i^v branching fraction from the likelihood fit, which depends on 
the form-factor model used to describe the signal distributions. Results for T^hy are given for each 
of the five form-factor models in Table |[ We use r^o = 1.56 it 0.04 ps 

The likelihood fit results are also used to measure the distribution in B p£v decay. 
The large b civ backgrounds in LOLEP limit this measurement to the HILEP lepton-energy 
region. To determine the distribution, we consider events in the vr^vr^ and 7r"'"7r~ modes with 
|M(7r7r) — M{p)\ < 0.15 GeV/c^ and \/S.E\ < 0.5 GeV. The background contributions in each 
interval are estimated using the results of the likelihood fit. We count the number of events in the 
background samples and subtract them from the data. Below we discuss the determination of q'^ 
for each event and describe the extraction of the distribution. 
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FIG. 7. The resolution for signal Monte Carlo events. The curve shows a fit to the Monte 
Carlo events (shaded histogram) using two Gaussians. The narrow (a = 0.4 GeV^/c^) component 
is due to events where the correct hadronic mode is reconstructed. The broader {a = 2.6 GeV^/c^) 
component is from crossfeed events. We account for this smearing in our measurement of AT in 
bins of g^. The full q"^ range for B pii' is Ag^ 21 GeV'^/c^. 

If the B meson momentum vector were known, could be determined unambiguously from 
the energy of the daughter hadron. For B — > piu, 

q^ = Ml + Ml - 2EbEp + 2pB • p^. (16) 

We know the magnitude of the B momentum, but not its direction, which introduces an uncertainty 
in q^. Using the beam-energy constraint, we compute the minimum and maximum allowed q^ 
values for each p + 1 combination. These correspond to the kinematic configurations where the B 
momentum direction is closest to, or furthest from, the p momentum direction. The midpoint of 
the allowed q^ range is our estimate of q^. The q^ resolution function is shown in Fig. ^. 

To find the true q^ distribution from the data, we must determine the efficiency of our analysis 
requirements as a function of q^, including the smearing effects shown in Fig. ^. We perform this 
measurement in three equal bins of q^. Monte Carlo simulation is used to compute the efficiency of 
events generated with a particular q^ to be reconstructed at a different value of (p' . This efficiency 
matrix is computed for each of the five form- factor models that we consider; model dependence arises 
primarily from the efficiency of the Ei^ > 2.3 GeV requirement. The reconstructed distribution 
of the HILEP data and this efficiency matrix are used to determine the true q'^ distribution. 



V. FIT RESULTS 

The signal and background yields extracted from the fit in HILEP and LOLEP are shown 
in Tables [V and ^ for the vrvr and ir^n^ir^ modes. The signal yields for the combined tt^tt^ 



and tt'^tt modes in HILEP, where our sensitivity to the signal is best, are 123 it 18 (direct) 
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HI 


HI 


uj HI 


p± LO 


pO LO 


w LO 


T(4S) yield 


198 


621 


460 


2249 


7298 


8552 


e+e" qq bkg. 


63± 11 


248 ± 22 


250 ± 22 


127 ± 16 


398 ± 28 


437 ± 29 


b civ bkg. 


39 ± 5 


90 ± 7 


52 ± 5 


-1 r\ At 1 A ^ 

1941 ± 47 


6419 ± 114 


7804 lb 129 


fake lepton bkg. 


4±2 


13 ± 7 


10 ± 5 


11 ± 7 


35 ±21 


44 ±22 


6 utv bkg. 


25± 10 


57± 23 


43± 17 


97± 24 


222 ± 57 


236 ± 60 


B — > -kIv bkg. 


9±3 


20 ±6 


4± 1 


15 ±5 


32± 10 


6± 2 


Direct sig. 


56 ±8 


67± 10 


21 ±3 


61 ±9 


79± 12 


31 ±5 


Xfeed sig. 


13 ±2 


80± 12 


32 ±5 


23 ±3 


100 ± 15 


52 ±8 



TABLE IV. Summary of data yields for the p , p , and u modes with lepton energies between 
2.3 and 2.7 GeV (HILEP) and between 2.0 and 2.3 GeV (LOLEP). Yields of the background 
contributions are insensitive to the form-factor model used to model B piu events; the yields 
presented in this table were obtained using the UKQCD model. The b — > uiiy background includes 
all B — > Xuii' modes except for p, uj, and vr. The crossfeed (Xfeed) signal contribution corresponds 
to events from one vector mode passing the selection cuts of another vector mode {p^iu p^lu 
or pii' <-> Loiv) and is constrained to the signal in the fit. All errors are statistical only. The errors 
on the direct and crossfeed signal yields are completely correlated, as are those on the yields of 
B ttIi/ events. Errors on the b — > uiv background are completely correlated among modes, but 
not between HILEP and LOLEP. 

and 93 ± 14 (crossfeed), for a total of 216 ± 32 signal events. The five b civ normalization 
scale factors are consistent with Monte Carlo predictions to better than 5%. The two scale factor 
ratios that determine the lepton-energy distribution for the b — > civ contribution are within one 
standard deviation of unity. The b uiv (modes other than B p£v, B uj£v, and B — > ttIv) 
normalizations in the three lepton-energy bins agree with each other to better than 20%, within 
one standard deviation, indicating that the ISGW2 cocktail of 6 — > uiv modes is adequate, at least 
to within our sensitivity. 

We show projections of the fit in both signal- and background-dominated regions of the AE 
vs. M(7r-7r) distributions. Figure ^ shows HILEP and LOLEP projections onto M{'inr) for \AE\ < 
0.5 GeV and onto AE for |M(7r7r) — Af(p)| < 0.15 GeV/c^. We observe a significant B piv signal 
in the M{t:tt) and AE distributions in HILEP. Figure ^ shows the same distributions for sidebands 
with lAi^l > 0.5 GeV and \M{'ktt) — M{p)\ > 0.15 GeV/c^, where we expect much less signal. 
In both signal- and background-dominated distributions, we observe good agreement between the 
data and the fit projections. In the projections shown, the signal component is modeled using the 
ISGW2 form-factor model, and the continuum contribution has been subtracted. 

The b civ background is quite small in the HILEP plots but is dominant in all of the LOLEP 
projections. The peak in the LOLEP vrvr modes at large M(7r7r) is due to S — > Div, with D — > Kir, 
where the K has been misinterpreted as a pion. LOLEP projections in Fig. ^ show that the shape 
of the b civ Monte Carlo distribution describes the data well in regions of the fit where the 
b — > uiv contributions are small. In Fig. ^, the shape of the M{Tnr) distribution in the LOLEP 
lepton-energy bin indicates that the b uiv contributions are necessary to describe properly the 
data below the b civ lepton-energy endpoint region. We do not show fit projections for the data 
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/9± HI 


/9° HI 


u HI 


/9± LO 


p° LO 


u LO 


Signal efF. (ISGW2) 


0.040 


0.091 


0.032 


0.039 


0.094 


0.040 


Signal yield (ISGW2) 


58 ±9 


68± 10 


21 ±3 


55 ±9 


70± 11 


26 ±4 


Xfccd yield (ISGW2) 


14 ±2 


80± 12 


34 ±5 


20 ±3 


89± 13 


48 ±7 


Signal efF. (LCSR) 


0.026 


0.061 


0.022 


0.030 


0.076 


0.034 


Signal yield (LCSR) 


56 ±8 


66± 10 


21 ±3 


64± 10 


84± 13 


33 ±5 


Xfeed yield (LCSR) 


14 ±2 


80± 11 


32 ±5 


24 ±3 


105 ± 15 


54 ±8 


Signal eff. (UKQCD) 


0.031 


0.070 


0.025 


0.033 


0.083 


0.036 


Signal yield (UKQCD) 


56 ±8 


67± 10 


21 ±3 


61 ±9 


79± 12 


31 ±5 


Xfeed yield (UKQCD) 


14 ±2 


80 ±12 


33 ±5 


23 ±3 


100 ± 15 


52 ±8 


Signal eff. (Wise/Ligeti+E791) 


0.035 


0.080 


0.028 


0.037 


0.094 


0.040 


Signal yield (Wise/Ligeti+E791) 


56 ±8 


66± 10 


21 ±3 


60 ±9 


79± 12 


30 ±4 


Xfeed yield (Wise/Ligeti+E791) 


14 ±2 


80± 12 


32 ±5 


24 ±3 


102 ± 15 


52 ±8 


Signal eff. (Beyer /Melikhov) 


0.029 


0.067 


0.024 


0.033 


0.083 


0.036 


Signal yield (Beyer /Melikhov) 


55 ±8 


66± 10 


21 ±3 


61 ± 10 


81 ± 12 


31 ±5 


Xfeed yield (Beyer/Melikliov) 


14 ±2 


80± 12 


32 ±5 


24 ±3 


104 ± 15 


53 ±8 



TABLE V. Summary of direct and crossfeed (Xfeed) signal yields for the p^, , and uj modes 
with lepton energies between 2.3 and 2.7 GeV (HILEP) and between 2.0 and 2.3 GeV (LOLEP). 
Signal efficiencies are normalized to the number of events in the full lepton-cncrgy range. The 
crossfeed contribution corresponds to events from one vector mode passing the selection cuts of 
another vector mode {p^iv p^lv or plv <-> uoli') and is constrained to the signal in the fit. All 
errors are statistical only. Errors on the direct and crossfeed signal yields are completely correlated. 
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FIG. 8. Projections of the fit for the combined and p'^ modes for HILEP {E^ > 2.3 GeV, 
upper plots) and LOLEP (2.0 < < 2.3 GeV, lower plots): (a) M{7nr) for HILEP after a 
|A^| < 0.5 GeV cut; (b) AE for HILEP after a \M{TnT) - M{p)\ < 0.15 GeV/c^ cut; (c) M(7r7r) for 
LOLEP after a |A^| < 0.5 GeV cut; (d) AE for LOLEP after a \M{7nr)-M{p)\ < 0.15 GeV/c^ cut. 
In each plot, the points with error bars show the on-resonance data after continuum-background 
subtraction, while the histogram shows the projection of the fit. The contributions to the fit are 
the direct and crossfeed components of the signal (unshaded regions, above and below the dashed 
line, respectively); the background from b — > uii' non-signal modes (darkly-shaded region); and 
the background from b — > civ (lightly-shaded region). The b uiv background includes B — > vr^i^ 
contributions. 
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FIG. 9. Projections of the fit for the combined and modes for HILEP (Ec > 2.3 GeV, 
upper plots) and LOLEP (2.0 < Ee < 2.3 GeV, lower plots): (a) M(7r7r) for HILEP after a 
\AE\ > 0.5 GeV cut; (b) for HILEP after a \M{7nT) - M{p)\ > 0.15 GeV/c^ cut; (c) Mirnr) for 
LOLEP after a |A£;| > 0.5 GeV cut; (d) AE for LOLEP after a |M(7r7r)-M(,9)| > 0.15 GeV/c^ cut. 
In each plot, the points with error bars show the on-resonance data after continuum-background 
subtraction, while the histogram shows the projection of the fit. The contributions to the fit are 
the direct and crossfeed components of the signal (unshaded regions, above and below the dashed 
line, respectively); the background from b uiv non-signal modes (darkly-shaded region); and 
the background from h clu (lightly-shaded region). The h ulv background includes B — > vrfz^ 
contributions. 
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FIG. 10. Projections of the kinematic variables: (a) and (b) E^. Both distribu- 
tions are shown for the combined and modes for events with |Ai?| < 0.5 GeV and 
|M(7r7r) - M{p)\ < 0.15 GeV/c^. The distribution is shown only for events with > 2.3 GeV. 
In each plot, the points with error bars show the on-resonance data after continuum-background 
subtraction, while the histogram shows the projection of the fit. The contributions to the fit are 
the direct and crossfeed components of the signal (unshaded regions, above and below the dashed 
line, respectively); the background from b — > uiu non-signal modes (darkly-shaded region); and 
the background from b — > civ (lightly-shaded region). The b uiu background includes B — > ttIv 
contributions. 



with Ei < 2.0 GeV, but the fit agrees well with the data in this region. The b — > cii/ Monte Carlo 
simulates M(7r7r) and AE projections of the data well in both shape and normalization. 

Figure ^ shows the reconstructed and lepton-energy distributions for the vrvr modes with 
|Aii^| < 0.5 GeV and |M(7r7r) — M{p)\ < 0.15 GeV/c^ requirements. We observe good agreement 
between the data and the distribution predicted in HILEP by the form-factor models. The 
lepton-energy spectrum predicted from the fit results shows good agreement with the data over 
the full range of lepton-energy used {Eg > 1.7 GeV). A large excess in the lepton-energy spectrum, 
consistent with the fit projection, is observed above the continuum contribution in the region beyond 
the b — > ciu lepton energy endpoint. 

In the above projections, the vr^vr'^ and vr+vr^ modes have been combined. Figure 11 shows 
the M^TTir) distribution with a \AE\ < 0.5 GeV requirement for these modes individually where 
the ratio of r(S° p-i+i^) to T{B+ p^l+v) has been allowed to fioat in the fit. We find 
T{B^ p-i+u)/T{B+ p^l+u) = l.ltll, where the error is statistical only. This value is in 
good agreement with the isospin relation in Eq. 13, which is used to determine our final results. 

In Fig. 12a we show the M(7r+7r~7r'^) plot for the HILEP lepton-energy bin. We do not observe 
a significant B — > ooiv signal, but the fit describes the data well. 

Figure |l^ also shows the AE distributions for the HILEP (|T^) and LOLEP (p^d) combined 
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FIG. 11. Projections of the fit onto Af(7r7r) for HILEP {Eg > 2.3 GeV) with a \AE\ < 0.5 GeV 
requirement: (a) for the tt^tt^ mode and (b) for the vr^vr^ mode. In each plot, the points with 
error bars show the on-resonance data after continuum-background subtraction, while the histogram 
shows the projection of the fit. The contributions to the fit are the direct and crossfeed components 
of the signal (unshaded regions, above and below the dashed line, respectively); the background from 
b — > uiiy non-signal modes (darkly-shaded region); and the background from b civ (lightly-shaded 
region). The b — > uiv background includes B — > iriv contributions. 
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FF model 


AF ng-i) 


AF (/10-2 ns-i) 


AF (/10-2 ns-^) 




< g2 < 7 GeVVc^ 


7 < g2 < 14 GeV2/c4 


14 < g2 < 21 GeVVc^ 




Ee > 2.3 GeV 


> 2.3 GeV 


Ee > 2.3 GeV 


ISGW2 


1.2 ± 0.5 


1.4 ± 0.8 


2.7 ± 0.8 


LCSR 


1.3 ± 0.5 


1.4 ± 0.8 


2.7 ± 0.8 


UKQCD 


1.2 ± 0.5 


1.4 ± 0.8 


2.7 ± 0.8 


Wise/Ligeti+E791 


1.1 ± 0.5 


1.4 ± 0.8 


2.8 ± 0.8 


Beyer /Melikhov 


1.1 ± 0.5 


1.4 ± 0.8 


2.8 ± 0.8 



TABLE VI. Results for AT{B^ p-~i+v) in bins of for events with E^ > 2.3 GeV for each 
form-factor (FF) model. The errors are statistical only. 



vr^ and vr'^ modes and the cos 6e (|T2|c) distribution for the HILEP and LOLEP vr^ and vr'^ modes. 
Independent of model, 9^ is expected to have a sin'^ 9 distribution for B ttIv events. The 
HILEP B Triu modes are dominated hy B ^ piv backgrounds. We find B{B^ vr^^+z/) = 
(1.3 ± 0.4) X lO"'', where the quoted error is statistical only, consistent with the previous CLEO 
result 1]^. We do not quote a full B -kIv result as it is very sensitive to systematics related to 
the large B — > p^v backgrounds. 

Table VI shows the results of the AF measurement before the correction for the lepton-energy 
cut is performed. The small spread seen among form-factor models is due to differences in the 
background subtraction and the smearing correction. Results for the full lepton-energy range are 
computed using the predictions in Table || and are discussed in Sec. VII. 



VI. SYSTEMATIC ERRORS 



Table VII summarizes our systematic errors. The largest systematic errors are due to un- 
certainties on the h — > ciu and b uiv backgrounds, the dependence of the efficiencies of our 
selection criteria on the Monte Carlo modeling of the detector, and the possible contamination due 
to — > TTTi^v nonresonant events. We consider these dominant uncertainties in more detail here. 

The leptons in 5 ^ civ background events are primarily from B — > Div, B D*ii', and 
B J/ipX decays. The B — > D**l!.u contribution is small but is nevertheless included in the 
fit. For events with |A£;| < 0.5 GeV and |M(7r7r) - M{p)\ < 0.15 GeV/c^ for the tttt modes and 
|M(7r+7r~7r°) - 0.782| < 0.05 GeV/c^ for the 7r+7r"7r° mode, the Monte Carlo simulation predicts 
the following mix of 6 ^ civ modes in HILEP (LOLEP): 23% (17%) for B Div, 48% (76%) for 
B D*£v, 23% (1%) for B J/tl^X, and 5% (5%) for B D**iv. 

The b civ systematic error reflects the fit sensitivity to the shape of the AE vs. M(7r7r) 
distribution determined from the b civ Monte Carlo simulation. (Recall that the fitting method 
already allows the data to determine the overall normalization of the b — > civ background as well 
as the shape of its lepton-energy distribution.) To evaluate this systematic error, we vary the 
relative sizes of the background sources. For B Div and B — > D*iv variations are ±20%, well 
beyond the current experimental uncertainties in their branching fractions. For B D**£v and 
B J/tpX we make variations of ±40%. These uncertainties produce only a small effect in the 
signal branching fraction, about 2%. 



26 



CM 
O 

0) 



CO 



(0 

c 
> 

LU 



10 



I I I I I I I I I I I I I I I I 




(a) HILEP 







' ' ' ' 



> 
0) 



o 

o 20 

CM ^" 



(0 

c 



3280399-012 

I I I I I I I I I I I I I I I I I I I 



(b) HILEP 



0.5 1.0 1.5 

MCtttttt) (GeV/c^) 



2.0 




-2-10 1 
AE(GeV) 



100 



tn 
+-> 

c 
> 

LU 



I I I I I I I I I I I I I I I I I I I 

(c) HILEP + LOLEP 




-1.0 -0.5 



cos 0, 



-2-10 1 
AE(GeV) 



FIG. 12. Projections of the fit in the vr^vr vr*^ and combined vr^ and vr*^ modes in HILEP and 
LOLEP: (a) the HILEP M(7r+7r~7r°) distribution in the 7r+7r-7r° mode with |AS| < 0.5 GeV; (b) 
the HILEP IS.E distribution for the vr modes; (c) the cos 6i distribution for the combined HILEP 
and LOLEP tt modes; (d) the LOLEP AE distribution for the vr modes. Oi is the angle between the 
lepton momentum direction in the W rest frame and the W momentum direction in the B rest frame. 
In each plot, the points with error bars show the on-resonance data after continuum-background 
subtraction, while the histogram shows the projection of the fit. The contributions to the fit are 
the direct and crossfeed components of the signal (unshaded regions, above and below the dashed 
line, respectively); the background from b — > uii' non-signal modes (darkly-shaded region); and 
the background from b ciu (lightly-shaded region). For the B uiiy channels, the b — > uiv 
background includes B Triv contributions. For the B niv channels, the b uiu background 
includes B — > p{uj)li' contributions. 
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Svst.PTTiatir rontriHntion 


6B IB 


'-'1 ^ ub\l 1 ^ ub\ 


bVjV 


jr/r 


6r/r 


y 1 cxii^c \ vjc- V / o J 












kjlliiUlclljlUll Ui LltJljCCLUl 


+Q% 


mo /o 


my /o 


mu /o 


□10 /o 


U ^ K^K-iy IjUiil UlJol LlUll 


+2% 


+ 1% 


+2% 

m^ /O 


□^u /o 


mo /o 


U u.-ct' OUliipUolLlUll 


□lU /o 


mo /o 


m^ /O 


□ly /o 


+7% 
m ( /o 


Tti t,p EfTri-t pd 1 n m i n osi t v 


±2% 


±1% 


±2% 


±2% 


±2% 


Lepton identification 


±2% 


±1% 


±2% 


±2% 


±2% 


Fake lepton rate 


±1% 


±1% 


±1% 


±3% 


±2% 


Fit technique 


±5% 


±3% 


±5% 


±5% 


±5% 


/+-//00 


±2% 


±1% 


±3% 


±3% 


±3% 


i3 tt'kIv nonresonant 


-8% 


-4% 


-10% 


-10% 


-10% 


TB 


±1% 


±2% 


±3% 


±3% 


±3% 


Total systematic error 


+13% 
-15% 


+7% 
-8% 


+12% 

-16% 


+14% 
-17% 


+13% 
-16% 



(excluding model dep.) 



TABLE VII. Summary of systematic errors. 

The same method is used to evaluate the uncertainty due to the h — > uif simulation. (The 
fit method allows the data to determine the normalization of this background component in each 
lepton-energy bin.) We vary the relative contribution of each mode predicted by the ISGW2 model 
by ±50%. The systematic error assigned is the sum in quadrature of these variations. 

We also examine the sensitivity of the fit to the B nil' decay distributions. In contrast to 
b — > uiiy, where the normalization is allowed to vary in each lepton-energy bin of the fit, the shape 
of the B — > iriu lepton-energy distribution is fixed by the Monte Carlo prediction. We therefore 
vary the lepton-energy spectrum of the B ttIv Monte Carlo and use form factors for B Triv 
determined using several different methods (much as we have done for the B — > pii' component 
of the fit). In both cases, the B piv results are robust against changes to the B — > irii/ Monte 
Carlo. We include this uncertainty in the systematic error due to the b — > uii^ simulation. 

The computation of |pmiss| relies on the Monte Carlo simulation to adequately simulate the 
detector response for all charged tracks and clusters in the event, not just those used to form the 
p + i candidate. To evaluate our sensitivity to the details of our detector Monte Carlo simulation, 
we vary the input parameters of the Monte Carlo. These variations include conservative variations 
of the tracking efficiency, charged-track momentum resolution, Csl cluster identification efficiency, 
cluster energy resolution, and the simulation of the detector endcaps. In addition, we examine our 
sensitivity to the number of mesons produced and the accompanying detector response. 

Finally, we assign a systematic error associated with a nonresonant B milu contribution. 
Like signal events, nonresonant events would have a IS.E distribution centered around zero. The tttt 
invariant-mass distribution, however, would be somewhat different. Although we do not know how 
to describe the M(7r7r) distribution of such a contribution, other properties help us to distinguish 
resonant B piv from nonresonant B — > nniv. The isospin (/) of the hadronic system in a 
B Tinlv decay must be either / = or / = 1. For / = 1, where the relative production of 
TT^vr'^rvr^Tr^iTr'^vr'^ is 2:1:0, the relative orbital angular momentum (L) must be odd. The L = 1 
contribution is dominated by the p resonance [||. Contributions from L = 3, 5, .... are suppressed. 
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For / = 0, the relative production of 7r^7r'^:7r"*"7r~:7r'^7r'' is 0:2:1, distinct from / = 1. An 7 = 
contribution will consist primarily of L = 0. In addition, the vr'^Tr'^ mode, which has no resonant 
contribution, is useful for constraining the size of any nonresonant B ttttIv contribution. 

Nonresonant events should also differ from B p£u events in their distribution. We expect 
nonresonant events to occur mainly at low g^, where the daughter u quark has a large momentum 
relative to the spectator quark. As shown in Table ||, the efficiency of the E£ > 2.3 GeV requirement 
is highest at large q'^. Thus, we expect to preferentially select resonant B piv events. 

To assign a systematic error, we include a possible nonresonant contribution in our fit, using 
two sets of assumptions for the M(7r7r) and distribution. First, we consider a M(7r7r) distribution 
with a broad Breit-Wigner shape having M = 0.8 GeV/c^ and F = 1.0 GeV/c^. To simulate a 
q^ distribution that is peaked at low g^, we use B — > irlu form factors from the ISGW2 model. 
The second set of parameters is designed such that the nonresonant contribution is very similar to 
resonant B pii/ events. We use B — > piu form factors (again from ISGW2) and a p Breit-Wigner 
for the parent M('inr) distribution. 

Using these two sets of assumptions, we repeat our likelihood fit with the additional freedom 
of a possible nonresonant B — > TTTriv contribution. As described above, we assume that this 
contribution is primarily / = 0. A 20% 1 = 1 contribution accounts for a possible L = 3 component. 
Additionally, we perform the fit with and without the tt^tt^ mode. 

In all fits, we find a nonresonant component consistent with zero. The systematic error is 
assigned from the increase in the statistical uncertainty on the B — > piv yield due to the correlation 
with the nonresonant fit component. This systematic error is one-sided, as our nominal fit assumes 
no nonresonant contribution. The associated uncertainty is 8%, which we regard as conservative. 



VII. B{B plv), \Vub\, AND AF RESULTS 

To extract the B — > piv branching fraction from our measured yields, we must extrapolate to 
the full lepton-energy range. We use form-factor models, as described in Section ||, to determine 
the signal efficiency. The fitted yields themselves depend only slightly on the set of form factors 
used. The signal efficiency and yield for each form-factor model are presented in Table 0. In 
addition, to extract \Vuh\-, we use the value of Fthy for each form-factor model, along with tb, to 



relate the branching fraction result to \ Vub\, as given by Eq. 15. Results for B{B ^ piiy) and \Vub\ 
are presented in Fig. |l^. 

The final values for B{B piu) and \Vub\ are the averages of the results obtained using the 
five form-factor models. We assign a systematic error to account for the substantial spread in 
results among the form- factor models. For the measurement of B{B ^ piu) this error is assigned 
to be 1/2 the full spread among the five form-factor model results. This uncertainty reflects our 
sensitivity to the different shapes of kinematic distributions predicted by the models, (mostly due 
to the Ei > 2.3 GeV acceptance. For \Vub\, we must also include an error due to the uncertainty 
on the computation of Fthy, relecting the different normalizations of the form-factor models. The 
models quote errors on Fthy between 17% and 50%. We have therefore assigned a 30% error on 
Tthy (corresponding to a 15% error on rather than 15%, which is the spread in predictions 

among models. We find 

B{B^ p-tv) = (2.69 ± 0.4lt°j5 ± 0.50) x 10-^ 

\Vub\ = (3.23 ± 0.24^0:1 ± 0.58) x 10~^ (17) 
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2.18±0.34*8;i§ 
3.19±0.48!§j^ 
2.79 ± 0.42*8 j| 
2.43 ±0.37^81 
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2.69 ±0.41 ±0.50 
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FIG. 13. B{B^ — > p~i~^u) and \Vub\ results. The errors are statistical, systematic, and theoret- 
ical (on the averages), successively combined in quadrature. For the branching fraction measure- 
ment, the theoretical error is taken to be 1/2 of the full spread of results. For the \ Vub\ measurement 
there is an additional contribution to the theoretical uncertainty due to the determination of Fthy 

The errors are statistical, systematic, and theoretical. The dominant uncertainty on \Vub\ arises 
from the theoretical error on the normalization, Fthy This 30% uncertainty (corresponding to a 
15% error on \Vub\) is independent of the method used to measure B{B pii^), and it is larger 
than the statistical error on \Vub\ (6%) or the model dependence of the detection efficiency (19% 
on e or 9% on 

Our results for AF in bins of are shown in Fig. 14. As for the branching fraction and \ Vub\ 



results, we compute an average over form-factor models and assign the theoretical uncertainty to be 
one-half the full spread in results. For these measurements, the model dependence comes primarily 
from the variation in the efficiency of the HILEP lepton-energy requirement. We find 



AF(0 < g2 < 7 GeVVc^) 
Ar(7 < g2 < 14 GeVVc^) 
Ar(14 < g2 < 21 GeVVc^) 



(7.6 ± S.Otii ± 3.0) X 10^2 ns-^ 
(4.8 ± 2.9+[5:| ± 0.7) X 10^2 ns-^ 



(7.1 ± 2.ll?;? ±0.6) X 10" 



ns 



(18) 



where the errors are statistical, systematic, and theoretical. Because the form-factor models predict 
nearly the same q'^ distribution at large lepton energy, we cannot distinguish among them, as shown 
in Fig. |l^. The models do, however, agree well with the q^ distribution seen in the data. At high 
q^, our lepton-energy requirement covers a large fraction of the allowed lepton-energy range (as 
shown in Fig. ^), and we are able to measure a partial rate with a relatively small theoretical error. 

Finally, we have formed an average for B{B plv) and \Vub\ with the previously published 
CLEO exclusive h — > ulv result |jl]. The two methods share only a small fraction of events and 
are therefore essentially statistically independent. The published result has been updated to con- 
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< q2 < 7 GeV^/c"* 
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FIG. 14. The partial width (AT) in bins of q^. The errors are statistical, systematic, and 
theoretical (on the averages), successively combined in quadrature. The theoretical error on each 
measurement is taken to be 1/2 of the full spread of results for the form-factor models. 
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FIG. 15. Comparison of measured AF distribution (points with error bars) with expectations 
from the form-factor models (curves) (a) for the Ei > 2.3 GeV region and (b) after the data have 
been extrapolated over the full Ei range. The errors on the points are statistical, systematic, and 
theoretical, successively combined in quadrature. Because the form-factor models predict nearly 
the same q'^ distribution at large lepton energy (a), we cannot distinguish between them. 
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FIG. 16. Branching fraction and \ Vub\ results averaged with the previous CLEO analysis. The 
errors are statistical, systematic, and theoretical (on the averages), successively combined in quadra- 
ture. For the branching fraction measurement, the theoretical error is taken to be 1/2 of the full 
spread of results. For the \ Vub\ measurement there is an additional contribution to the theoretical 
uncertainty due to the determination of Fthy ■ 



sider the same set of form-factor models as used in this paper. (These results are described in 
Appendix ^.) A weight has been assigned to each analysis such as to minimize the total error 
(statistical, systematic, and theoretical) on the average. Because Ref. |[l| extracts using both 
B — > piu and B — > niu results, we perform the average separately for the branching fraction and 
\Vub\ measurements. Figure ^ shows the resulting averages accounting for correlated systematic 
errors between the two methods. Averaging these results over form-factor models, as described 
above, we obtain 

B{B^ p-l+v) = (2.57 ± 0.29i[!:^^ ± 0.41) x 10"^, 

\Vub\ = (3.25 ± 0.14t[}:^^ ± 0.55) x 10"^ (19) 

where the errors are statistical, systematic, and theoretical. 



VIII. CONCLUSIONS AND OUTLOOK 

We have performed a measurement of B{B —>■ piu), \Vub\, and the distribution in — > piu 
decay using a data sample of approximately 3.3 x 10^ BB pairs. Using leptons near the b — > ciu 
lepton-energy endpoint, we find 

B{B^ p-tu) = (2.69 ± 0.4lt°;^5 ± o.50) x 10~^ 

\Vub\ = (3.23 ± 0.24^0:1 ± 0.58) x 10"^ (20) 
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where the errors are statistical, systematic, and theoretical. The B{B plv) result confirms the 
previous CLEO measurement and has a comparable statistical precision. This result is statisti- 
cally independent from the previous CLEO result and has a somewhat smaller systematic error. 
Averaging the measurements, we find 

B{B^ p-£+iy) = (2.57 ± 0.29l[|j^ ± 0.41) x 10-^ 

\Vub\ = (3.25 ± O.Utoil ± 0.55) x IQ-^. (21) 

These values represent the current best CLEO results based on exclusive b ulu measurements. 
For the branching fraction, the experimental a^id theoretical uncertainties (±16%) are 

comparable. For \Vub\- however, the experimental uncertainties (I^q^) substantially smaller 
than the estimated theoretical error (±17%). The theoretical error on \Vu,b\ contains a contribution 
from the uncertainty on the detection efficiency (±8%) and one from the overall normalization 
(±15%). 

We have also measured the distribution in S — piu decay in three bins. We find 

Ar(0 <q^ <7 GeVVc^) = (7.6 ± S.Ot?! ± 3.0) x 10"^ ns'^, 
Ar(7 < g2 < 14 GeVVc^) = (4.8 ± 2.9toi ± 0-7) x lO'^ ns'^, 
Ar(14 < < 21 GeVVc^) = (7.1 ± 2.1;?;? ± 0.6) x 10"^ ns'^. (22) 

Because we are sensitive to a large fraction of the allowed lepton-energy region at high g^, we 
measure the partial rate for 14 < < 21 GeV^/c^ with a relatively small theoretical uncertainty. 
This result is promising for future analyses oi B ^ p£u that use the high lepton-energy region to 
determine \Vub\, if theoretical predictions for AFthy at high can be made with good precision. 
For high lepton energies, the predicted shapes of the q"^ distributions arc virtually identical for 
the different form-factor models. The models, therefore, cannot be distinguished on the basis of 
measurements in this lepton-energy region alone. 

To determine |V^b| more precisely using the full q^ range, it is important to include leptons whose 
energy is below the lepton-energy endpoint region, so that the full distribution of B ^ p£u decay 
can be measured. Experimental results on the g^ distribution for leptons with E( < 2.0 GeV would 
help to improve form- factor models, and thus the measurement of \Vub\- 
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APPENDIX A: UPDATE OF PREVIOUS CLEO RESULTS 



In this Appendix we present updated results of the previously published CLEO analysis of 
B piu and B tt£u, which used 2.84 x lO^BB pairs M. These results are used in Section [VI I| 
to compute average values of B{B ^ piu) and \Vub\. 
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piu FF model 



niu FF model 



B{plv) (/10~4) B{t:Iv) (/lO- 



2.2 ±0.4l^:^ 
2.8 ±0.5+[j| 

2.5 ±o.5;[|:f 

2.2 ±0.4t[}-^ 
-0:5 



IK 



ub\ 



ISGW2 
LCSR 
UKQCD 
Wise/Ligeti+E791 
Beyer /Melikhov 



ISGW2 
QCD SR 
QCD SR 
QCD SR 
Beyer /Melikhov 



2.5 ±0.5] 



0.7 



2.0 ± 0.5 ±0.3 

1.6 ± 0.4 ±0.3 

1.7 ± 0.4 ±0.3 

1.8 ± 0.4 ±0.3 
1.6 ± 0.4 ±0.3 



3.3 ±0.2; 

3.4 ±0.2: 
3.3 ±0.2: 
3.0 ±0.2: 
3.3 ±0.2l:[Jj 



-0.4 
-0.3 
-0.4 
-0.3 
-0.4 
-0.3 
-0.4 



TABLE VIII. Updated results of Ref. § for B{B^ p'^^i^), B{B^ -k-I+u), and \Vuh\. The 
errors are statistical and systematic. 



We have updated these results to use form-factor models for B piu and B — > Triv decay 
developed since Ref. [||. To compute B{B — > p£u) and \Vub\, a form-factor model is needed to 
describe both B piu and B — > irlv decay. In general, we use predictions from the same 
authors to describe both sets of form factors when both are predicted [^,10|. For the fits using 
the UKQCD, LCSR, and Wise/Ligeti±E791 predictions to describe B p£u events, we use the 
B ttIv predictions from Khodjamirian et al. [|l8| using QCD sum rules (QCD SR). 

The fit results have been updated to use these form-factor models as well as the tb results 
used in this paper. No other changes from Ref. Q have been made. Results for B{B — > piu), 
B{B Triu), and \Vub\ are shown in Table VIII . 
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